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Abstract 

One-parameter family of discrete-time quantum-walk models on the square lattice, which includes 
the Grover-walk model as a special case, is analytically studied. Convergence in the long-time limit 
t — ► oo of all joint moments of two components of walker's pseudovelocity, Xt/t and Yf/t, is proved 
and the probability density of limit distribution is derived. Dependence of the two-dimensional limit 
density function on the parameter of quantum coin and initial four-component qudit of quantum 
walker is determined. Symmetry of limit distribution on a plane and localization around the origin 
are completely controlled. Comparison with numerical results of direct computer-simulations is 
also shown. 
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I. INTRODUCTION 



Quantum walks are expected to provide mathematical^irmdels^for quantum algorithms, 
which could be used in quantum computers in the future^ 
atic study of quantization of random walks is not old 
have been well studied and mathematical properties are clarified 




, |5| . Though the system- 
9j, one-dimensional models 



10 



For example, 



convergence of all moments of pseudovelocity in the long-time limit was proved for the 
standard two-component quantum-walk model and the weak limit-theorem is established 
13, Q, [ijj]. The weak limit-theorem was generalized for the multi-component quantum- 
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walk models associated with rotation matrices [16l . [1 

One of t re recent topics of quantum walks is systematic study of higher dimensional 



models [14], \m, 



19 



20 



21 



22| . Among them the Grover-wa 



23 



k model has been extensively 



2% 123, M, Ui- Inui et aL[28| 



studied, since it is related to Grover's search algorithm 
studied the two-dimensional Grover-walk model analytically and clarified an interesting 
phenomenon called localization 29(. In two dimensions effect of random environment on 
quantum systems is non-trivial and decoherence in two-dimensional quantum walks generated 
by broken-line-type noise was studied by Oliveira et al. 
We noted that at the end of the paper by Inui et al 



3Q|. 



a one-parameter family of two- 
dimensional quantum-walk models was introduced, which includes the Grover walk as a 
special case; with the parameter p = 1/2 of a quantum coin. In general the quantum walker 
on the square lattice, which hops to one of the four nearest-neighbor sites at each time step, 
is described by a four-component wave function. In the present paper, we will determine 
the dependence of long-time behavior of quantum walker both on the parameter p and a 
four- component initial wave function (four-component qudit) completely and establish the 
weak limit-theorem for the family of two-dimensional models. 

This paper is organized as follows. In Sec. II we define the discrete-time two-dimensional 
quantum-walk models. By calculating the eigenvalues and eigenvectors of the time-evolution 
matrix of quantum walk in the wave-number space, long-time behavior of joint moments of 
x and y components of pseudovelocity is analyzed in Sec. III. There the weak limit-theorem 
for the two-dimensional models is proved and dependence of the limit distributions of pseu- 
dovelocities on the parameter p of quantum coin and on an initial qudit of walker is clarified. 
In order to demonstrate the usefulness of our results to control the long-time behavior of 



2 



quantum walks, we show pairs of figures of direct computer-simulation results and of ob- 
tained limit distributions in Sec. IV. Using the results we can discuss symmetry of limit 
distributions on a plane systematically depending on the parameter p and initial qudits of 
walker. Concluding remarks are given in Sec.V. Appendix A is used to show calculation of 
some integrals. 



II. TWO-DIMENSIONAL QUANTUM WALK MODELS 



A. General setting on the square lattice 



We begin with defining the two-dimensional discrete-time quantum walk on the square 
lattice Z 2 = {(x,y) : x,y G Z}, where Z denotes a set of all integers Z = 
{• • • , — 2, — 1, 0, 1, 2, • • •}. Corresponding to the fact that there are four nearest-neighbor 
sites for each site (x,y) G Z 2 , we assign a four-component wave function 



*(x,y,t) 



^ 4>i(x,y,t) ^ 

^s(x,y,t) 
yip 4 (x,y,t) J 

to a quantum walker, each component of which is a complex function of location (x, y) G Z 2 
and discrete time t — 0, 1, 2, • • •. A quantum coin will be given by a 4 x 4 unitary matrix, 
A = (Ajk)j k=1 , and a spatial shift-operator on Z 2 is represented in the wave-number space 
(k x , k y ) G [— 7r, 7r) 2 by a matrix 



S(k x , k y ) — 






V 





















/ 



where i = We assume that at the initial time t — the walker is located at the origin 

with a four-component qudit T o — (<?i? Q2, 93, Q4) G C 4 , Yfj=\ \lj\ 2 — 1- 111 the present paper, 
the transpose of vector/matrix is denoted by putting a superscript T on the left, and R and 
C denote the sets of all real and complex numbers, respectively. Let 



V(k x , ky) — S(k x , ky)A. 



(1) 



Then, in the wave-number space, the wave function of the walker at time t is given by 

if(k x , ky, t) = (V(k x , kyj) t = 0, 1, 2, ■ • • . 



(2) 



Time evolution in the real space Z 2 is then obtained by performing the Fourier transforma- 
tion 

J — 7r Z7T J-vr Z7T 

Note that the inverse Fourier transformation should be 



i>{k x ,k y ,t) = J2 ^(x,y,t)e-^ x+k y^. 

Now the stochastic process of two-dimensional quantum walk is defined on Z 2 as follows. 
Let X t and Y t be x and ^-coordinate of the position of the walker at time t, respectively. 
The probability that we find the walker at site (x, y) G Z 2 at time t is given by 

P(x, y, t) = Prob((X t , Y t ) = (x, y)) = &(x, y, t)*(x, y, t), (3) 

where ^(x, y, t) = T ^(x, y, t) is the hermitian conjugate of ^(x, y, t). The joint moments of 
X t and Yt are given by 



= x a y^P(x,y,t) 



._d_ 

3kn, 



^(t) t), 



(4) 



for a, (3 = 0,1,2, 



B. Generalized Grover walks 



Inui et al. 



281 ] introduced a one-parameter family of quantum-walk models on Z 2 as a 



generalization of Grover model by specifying the quantum coin as 



' -p q Jpq Jpq^ 



A 



q -p ^jpq ^jpq 
pq \Jpq~ —q p 
pq yfpq p -q J 



q=l-p, 



(5) 



where p G (0, 1). When p = 1/2, A is reduced to the quantum-coin matrix used to generate 
the Grover walk on Z 2 . In general the generator of the process ([1]) is given as 



/ 



V (k x , ky 



-pe 



I k'x 



qe 



qe 
-pe 



pqe lky 
\ \Jpqe~ lky 



pqe lky 



pqe 



q = 1 — p,0 < p < 1. 



pqe 



pqe 
-qe 



iky 



ik-r 



pqe 

pqe ~ lkx 
-qe~ tky 



\ 



III. LIMIT DISTRIBUTION IN t -» oo 



A. Calculation of moments and their long-time limits 



(6) 



In order to analyze the long-time behavior of the present two-dimensional quantum walks, 
we use the method originally given by Grimmett et al. 1141, which has been developed in 

nnn 

mm na. it is easy to diagonalize the time-evolution matrix ([6]). The four eigenvalues 
are obtained as 



\ x = 1, A 2 = -1, A 3 = e^'H A 4 = e 



■iu}(k x ,k y ) 



where w(k x , k y ) is determined by the equation 



cosu(k x , k y ) = — (p cos k x + q cos k y ) 



(7) 



The eigenvectors corresponding to the eigenvalues Xj, 1 < j < 4, are given by the following 
column vectors 



Vj(k x ,k y ) = Nj 



(8) 



' qie^Xj + l)(e ifc «Aj + l)(e _< *»A i + 1) ^ 
q(e iky \j + l)(e- ik *\j + l)(e"^A j + 1) 
Jpq(e iky \j + l)(e^ fe A j + l)(e ife A i + 1) 
I y/pq(e~ ik "*j + l)(e ik *\j + l)(e"^A j + 1) J 

with appropriate normalization factors Nj,l < j < 4. Define the 4x4 unitary matrix 
R(k x ,k y ) = (vi, v 2 , v 3 , v 4 ) from the four column vectors (jHJ). Then the time-evolution 
matrix (jSJ) is diagonalized, and by the unitarity of R(k x , k y ), R^(k x , k y ) = [^(A^, ky)]' 1 , (j2J) 



5 



is written as 



1 \\ ^ 



o xi o o 



A* 3 



A 



y)9o 



4/ 



(Aj) VjCj(k x , k y ), 



i=i 



where Cj(k x , k y ) = v](/c x , k y )(j)o. For a, /3 = 1, 2, • • -, we see 



^(k x ,k y ,t) 



dk 7 



+ 



^ du(k x ,k y ) J (M^k^ p ( a 4 )V 4 (^, yc 4 ^, fc,)^ + o^- 1 ; 



since Ai = 1 and A2 = —1 are independent of /c x , fc^. Since R(k x , k y ) is unitary, its column 
vectors make a set of orthonormal vectors; v^(A; x , k y )v m >(k x , k y ) = 5 mm >. Then we have 

= {(-ir + e\c 3 (k x , k y )\ 2 + \c 4 (k x , Ag| 2 } ( ^ y )^ a+/3 + o^- 1 ). 

The pseudovelocity of quantum walker at time t is defined as 



V 



t_ l T'T 



t = 1,2,3, 



(9) 



Eq.(J3]) gives the following expression for joint moments of x and y components of pseudove- 
locity, (X t /t) a (Y t /t) 13 , in the long-time limit 



lim 

t^co \ \ t 



X t \ a fY t ^ s 



t 



T (lki r ,iky {(-ir +0 \c 3 (k x ,k y )\ 2 +\c,(k x ,k y )\ 2 } 



x 



' duu(k x ,k y )\ a fdu(k x ,k y ) 







\ J V dky J 

Here from (JTj) we have u(k x , k y ) = arccos { — (pcos k x + qcosk y )} and then 

du)(k x ,k y ) psink x 



dk x 

du(k x , k y 
dk„ 



1 — (p cos k x + q cos ky) 2 

q sin ky 
1 — (p cos k x + q cos fcy) 2 



by the formula (d/dx) arccosx = =Fl/y 1 — x 2 . 

We change the variable of integral from k x , k y to v x , v y by 

p sin k x 

?;„ = — — 
1 — (p cos k x + q cos /Cj,) 2 

g sin k y 



1 — (p cos + q cos fcy) 



(10) 



It should be noted that this map (k x , k y ) € [— ir, n) 2 \— > (u„, is one-to-two and the image 



is a union of interior points of an ellipse 



v 2 v 2 
^ + ^<1 



(11) 
P Q 

and the four points {(p, q), (p, —q), (— p, q), (—p, — <?)}• We found that the following relations 
are derived from ffTUl) . 



sin k r 



cos k x 



sin L, 



cos fc,, 



^ x ^/pq-qv 2 -pv 2 



P\J(V~X + Vy + 1)(V X ~Vy + l)(V X +Vy~ 1)(V X -Vy-1) 

(l + q)v 2 x + pv 2 y -p 



Pj(V X +Vy+ 1)(V X ~Vy + 1){V X + Vy~ l)(V X -Vy-1) 



2v y Jpq - qv 2 x -pvl 



q\j (Vx + Vy + l)(V x ~Vy+ 1)(V X +Vy~ 1)(V X - Vy-1) 

1 v l + ~ 1 



(12) 



q\J (V X +Vy+ 1)(V X ~Vy + 1)(V X + Vy- l) {V X ~ Vy ~ 1 ) 

They are useful to calculate the Jacobian associated with the inverse map (v x , v y ) i— > yk x , ky) 



and we have obtained 



J 



dv x /dk x dv x /dk y 
dv y /dk x dvy/dky 

(V X + Vy + 1)(V X -Vy + 1)(V X +Vy - l)(V X ~ V y ~ l) 



If we assume that by this change of variable Cj(k x , k y ) are replaced by Cj(v x ,v y ), j = 3,4, 
the integral is written as 

.im,i^m w 



(13) 



i^oc \ \ t J \t 

1,00 dv, r f°° dv„ 1 



2tt 



OO ZjVI j —CO 
OO /'OO 



"—j {\Cs(v x ,v y )\ 2 + (-l) Q ^|(7 4 K,^)| 2 }^l M/p+ui/9<1} 



2tt J 



OO •/ — oo 



dv x / dv y v x v^p(v x ,Vy)M( 



where denotes the indicator function of a condition fi; = 1 if SI is satisfied and 
= otherwise. Here ^ p {v x ,v y ) is given by 

2 

^ Vy) = -n\v x + v y + l)(v x -v y + l){v x + v y - l)(v x -v y -l) ^.WKV' (14) 

since we can confirm that (v x + v y + l)(v x — v y + l)(v x + v y — l)(v x — v y — 1) > 0, when 
v x /p+v y / q < 1, q — l—p, < p < 1. This function fj, p (v x , v y ) gives the fundamental density- 
function for long-time limit distribution of pseudovelocity (see Appendix A). Figure [H shows 
it when p — 1/4. It should be noted that the fundamental density-function fi p (v x , v y ) depends 
on the parameter p but does not on an initial qudit T O — (<7i, ?2> Qs, Qi)- The dependence 
on an initial qudit is expressed by the weight function A4(v x , v y ) given below. 




FIG. 1: (Color online) The two-dimensional fundamental density-function )J, p (v x ,v y ) of limit dis- 
tribution of pseudovelocities, when p = 1/4. 
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B. Weight function Ai(v x ,v y ) 

Using (H2D, the weight function Ml(v x , v y ) is explicitly determined as follows; 

M(v x , v y ) = Mi + M 2 v x + M 3 v y + M A v 2 x + M b v 2 + M 6 v x v y 



(15) 



with 



Mi 
Mi 

M 3 



- + Re(gig 2 + g3<?4), 
-(M 2 - l^l 2 ) + -j=Re{qiq 3 + qm - q 2 q 3 ~ 42(74), 

p 



■(ifcl'-lfcp 



Re(<?i<?3 - q^ 4 + q 2 q 3 - q 2 q±), 



M 4 = \(\qi\ 2 + |<? 2 | 2 - N 2 - M 2 ) - — Re(M 2 ) - Refofc) 
2 v ' p 

Re(gig 3 + qiq A + q 2 q 3 + q 2 q A ), 

pq 

-l(\Qi\ 2 + N 2 - N 2 - M 2 ) - Re(gig 2 ) - i±*Re( ft &) 
L q 

p 



Mr 



Re(<?i<? 3 + q x q A + q 2 q 3 + q 2 q 4 ), 



M t 



pq 
1 



Re(gig 3 - q x q± - q 2 q 3 + q 2 q 4 ), 



(16) 



where Re(z) denotes the real part of z G C and z denotes the complex conjugate of z G C. 
The weight function defines the following real symmetric matrices M n , through the relations 

M r 



i„ = </)jM n o , 1 < n < 6, 



M, 



^ 1 1 0^ 
110 
11 
11 

^00 




^ 2\fpq —q—q^ 



M 9 



—2^pq q q 
-q q 
1 -g g 



2Jpq 



—p p 
—p p 
—p —p 2^pq 
y pp -2^Jpq 



-1 


111 


q 


q 


P 


Vm 


Vpq 


1+9 


-1 


q 


q 


P 




Vpq 


Vpq 


1 


q 


1 


1 




Vpq 






q 


q 


1 


1 


Vm 


Vpq 
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M, 



p p \ 

ffq yfpq 

P P 

I 1+2 



VpI sfm q 

p p 1+2 _i 
Vpq Vm i 



M fi 



2Jpq 



J 






-1 1 



-1 1 
1 -1 





V 



1-10 



/ 



Such matrix representations will be useful, when we generalize the present results to other 
models, whose quantum coins are given by larger matrices 171 ]. 

The integral /f^ dv x /f^ dv y fj, p (v x , v y )A4(v x , v y ) is generally less than one, since the con- 
tributions from the eigenvalues Ai and A2 have not been included. The difference 



A = 1 - 



dV x / dVyjJ p (V X ,Vy)M(V X ,Vy 



(17) 



gives the weight of a point mass at v x 



in the distribution. That is, A gives the 



probability of localization around the origin of the present two-dimensional quantum walks 



IS 



(see Sec. III. D below). 



C. Weak limit-theorem and symmetry of limit distribution 

The result is summarized as the following limit theorem. 
Theorem Let 

V{V X ,Vy) = flp(V X ,Vy)M(V X ,Vy) + A5 (V x ) 5 (Vy ) , (18) 

where fJ- p (v x ,v y ), A4(v x ,v y ), and A are given by (fUjl . ffT5|) with fTl6l) . and ffl7|) . respectively, 
and 8(z) denotes Dirac's delta function. Then 

hm l(^f] = r dv x r dv y vy v v(v x ,vy) (19) 



00 j — 00 



t^oo \ \ t J \t 

for all a, (3 = 0, 1, 2, 

As mentioned in an earlier paper [l6|, distribution of quantum walks itself does not 
converge in the long-time limit, since time evolution of quantum system is simply given by 
a unitary transformation. The above theorem is regarded as a weak limit-theorem in the 
sense that, if we evaluate moments of pseudovelocity in oscillatory distributions of realized 
quantum walks, the results shall be converge to the values calculated by the formula ( fl9l) 
with the density function ( fl8l) in t — > 00. If we integrate u(v x ,v y ) over any region Dona 
plane R 2 , then we obtain the probability that the pseudovelocity V t = (X t /t,Y t /t) G D in 
the t — * 00 limit. 
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The polynomial form of (fT5j) leads to the following classification of symmetry realized in 
the limit distribution. 

(i) When M3 = Mq = 0, the limit of probability density v{v x ,Vy) has the reflection 
symmetry for the t^-axis; v{y x , — v y ) = is(v x ,v y ). 

(ii) When M2 = M$ = 0, the limit of probability density u(v x ,v y ) has the reflection 
symmetry for the tyaxis; v(—v x ,v y ) = z/(t> x ,w y ). 

(hi) When M2 — M$ = Me = 0> the limit of probability density u(v x ,v y ) has the 
reflection symmetries both for the t^-axis and the tyaxis; u{v X) —v y ) = u(—v x ,v y ) = 
u(v x ,v y ). 

(iv) When M. 2 = M3 = 0, the limit of probability density u(v x ,v y ) has the bi- 
rotational symmetry for the t^-axis, which is perpendicular both to v x - and 

V(~V X , ~Vy) = l>(V X ,V y ). 



D. Localization probability around the origin 

By symmetry of the fundamental density-function ffT4l) . (fl7|) with ( fl5|) becomes 

A = 1 - Mi - M A K X - M 5 K y 

with 

/oo /*oo 
dv x / dv y [i v {v x ,v y )v\ 
-00 J — 00 



2 

XI 

00 J —00 
00 /-oo 



K y = / dv x I dvyHp(v X} v y )v y . 

As shown in Appendix A, these integrals are readily performed and we obtain the following 
explicit expression for the probability of localization around the origin, 

2 2 
A = 1 — Mi (arcsin Jp — Jpq)Mi (arcsin Jq — ^/pq)M^. (20) 

7T 7T 

The localization probability A is a function of the parameter p e (0, 1) and an initial 
four- component qudit T O = (<?i, Q2, Q3, Q4) £ C 4 , Y^=i \lj\ 2 — 1 through ([TBI) . For example, 
(EOD gives 

1 / ^ f 1 ^ 1 1 ] 

A = — (1 — 2 v /pg) < - arcsin y/p H — arcsin 



7T 



I P 9 a/P^ . 
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for r o = (1, 1, -1, -l)/2, and 

A = — (1 + 2Jpq) \ - arcsin Jp H — arcsin Jq 1 

tt [p q jpq) 

for T (j) = (1,1,1, l)/2, respectively, where q — 1 — p. As shown in Figl21 for T <p = 
(1, 1, — 1, — 1)/2, the localization probability A attains the minimum = for the Grover- 
walk model, p = q = 1/2, while for T 0o = (1,1,1,1) /2, it attains the maximum 
= 2(7r — 2)/ tt = 0.726 ■ • • for the Grover-walk model. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

P P 



FIG. 2: (Color online) Dependence of localization probability around the origin A on the parameter 
p G (0, 1). (a) The case T cf)o = (1, 1, -1, -l)/2. Whenp = 1/2 (the Grover-walk model), A = 0. (b) 
The case T (j) Q = (1, 1, 1, l)/2. When p = 1/2 (the Grover-walk model), A = 2(vr - 2)/vr = 0.726- • -. 



If we make the initial qudit depend on the parameter as 

0O =lV2'V2'-V2'-V2j' q = 1 ~ P ' (21) 

for example, then A = for M.\ = 1,M.4 = = 0, and thus the quantum walker is 
extended with probability one for all p £ (0, 1). 

It should be noted that A is defined as the intensity of Dirac's delta-function at the 
origin found in the limit density-function of pseudovelocity (see Eq. f|T8|) ) . It implies that 
A gives the probability that the quantum walker loses its velocity and stays around the 
starting point, i.e. the origin. Therefore, A is, in general, greater than the time-averaged 
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probability that the walker stays exactly at the starting point, P^, which was calculated 

Tj 



in [28] . For example, for the Grover-walk model with the initial qudit 0o = (1? 1> 1; l)/2, 



0.726' 



A = 2(vr - 2)/tt 



reported in Sec.V.C in [28]. 



as mentioned above, while P c 



2{(tt - 2)/tt} 2 = 0.264... as 



IV. COMPARISON WITH COMPUTER SIMULATIONS 

In order to demonstrate the validity of the above results, here we show comparison with 
numerical results of direct computer simulations [16]. In FigsJHJEl the left figures show the 
distribution of pseudovelocity V$ = (X t /t,Y t /t) at time step t = 30 numerically obtained 
by computer simulations and the right figures the long-time limits of probability densities 
v(v x , v y ) determined by our theorem. The four figures show the symmetries (i)-(iv) classified 
in Sec. III. C. In all of these four cases shown in Figsj3][6l A > and we can see a peak at 
the origin in each right figure (b), which indicates the contribution A5(v x )S(v y ) in the limit 
density-function ( TT8l . 




FIG. 3: (Color online) The case p = 1/4 and T cf) = (1, -1, 1, l)/2. Since M 3 = M 6 = in this 
case, the limit distribution has the reflection symmetry for the i? x -axis; v(v x ,—v y ) = u(v x ,v y ). 
(a) Distribution of pseudovelocity Vj = (X t /t,Y t /t) at time step t = 30 numerically obtained by 
computer simulation, (b) Probability density of limit distribution. 



We observe oscillatory behavior in distributions of V 4 = (X t /t,Y t /t) in computer simu- 
lations. In general, as the time step t increases, the frequency of oscillation becomes higher, 
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FIG. 4: (Color online) The case p = 1/4 and T 4> = (1, 1, 1, -l)/2. Since .M 2 = -M 6 = in this 
case, the limit distribution has the reflection symmetry for the v y -axis; v(—v x ,v y ) = v(v x ,v y ). 
(a) Distribution of pseudovelocity Vj = (Xt/t,Yt/t) at time step i = 30 numerically obtained by 
computer simulation, (b) Probability density of limit distribution. 




FIG. 5: (Color online) The case p = 1/4 and T 4> = (1, 1, 0, 0)/\/2. Since M 2 = M 3 = M 6 = 
in this case, the limit distribution has the reflection symmetries both for the v x -axis and the v y - 
axis; u(v x , -v y ) = v{ ~v x ,v y ) — u(y x ,Vy). (a) Distribution of pseudovelocity — (Xf/t,Yt/t) at 
time step t = 30 numerically obtained by computer simulation, (b) Probability density of limit 
distribution. 



but, if we smear out the oscillatory behavior, the averaged values of distribution shall be 
well- described by the density functions of limit distributions (I18p . which is the phenomenon 
implied by our weak limit-theorem 
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FIG. 6: (Color online) The case p = 1/4 and T <t> = (1, -1, -1, l)/2. Since M 2 = M 3 = 
in this case, the limit distribution has the bi-rotational symmetry for the v z -axis, which is per- 
pendicular both to v x - and i^-axes; v{—v X} —v y ) = v(v x ,Vy) (a) Distribution of pseudovelocity 
Vf = (Xf/t,Yt/t) at time step t = 30 numerically obtained by computer simulation, (b) Probabil- 
ity density of limit distribution. 



V. CONCLUDING REMARKS 

In general, quantum coins, which determine time-evolution of quantum walkers with 

n n 

spatial shift-operators, are given by unitary transformations [4J, |16j. The set of all N x N 
unitary matrices makes a group, the unitary group U(N), whose dimension is N 2 (see, 



for example, 31]). Though the determinant of unitary matrix is generally given by e t(p ,ip G 
[— 7r/2, vt/2), this global phase factor of quantum coin is irrelevant in calculating any moments 
of walker's positions in quantum-walk models [lj3|. For example, in the standard two- 
component (N = 2) quantum walks, the number of relevant parameters to specify a quantum 
coin is N 2 — 1 = 2 2 — 1 = 3 (Cayley-Klein parameters), and the dependence of limit 



distributions of pseudovelocities on the three parameters was completely determined 
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3, 



13l . Il5l . 1 1 €31 ] . In the present paper we have considered a one-parameter family of unitary 
matrices (jSJ) in U(4) as quantum coins. The present study should be extended to more 
general models, whose U(4)-quantum coins are fully controlled by 4 2 — 1 = 15 parameters. 

One of the motivations to study the present family of models in this paper is the fact 
that it contains the Grover walk on the plane. It will be interesting and important to 
derive limit distributions of pseudovelocities of quantum walkers on varietyof plane lattices 
different from the square lattices and in the higher- dimensional lattices [20J. For example, 
the quantum coin of the Grover walk in the D-dimensional hyper-cubic lattice is given by 



15 



the 2D x 2D orthogonal matrix = (A^) with the elements 



1/D-l, ifj = k 
1/D, iij^k. 



It is also an interesting problem to relate the present results to solutions of the continuous- 
time quantum- walk models on two-dimensional lattices j^]. 

At the end of the present paper, we refer to the fact that recent papers propose im- 



plementations of not on 



optical equipments [32 



cal lattices 



y 



one-dimensional but also two-dimensional quantum walks using 



35 



33], ion-trap systems [3J], and ultra-cold Rydberg atoms in opti- 
36| . We hope that combinations of experiments and theoretical works of 



quantum physics will make significant contribution to development of quantum informatics. 
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APPENDIX A: ON INTEGRALS 

Consider the integral 

/OO /"OO j 
dv x / dv y l {v 2j +v 2, <1} - ■ — — — — ■ — — 
-OO J-OO (V X + Vy + 1){V X - Vy + 1)(V X + Vy ~ 1){V X ~ Vy ~ 1) 

with p + q = 1, p, q > 0. Let 

(ad 

Then 

Mr) 



I = — 2 4 z v /pg / dr 
Jo 

with a contour integral on a complex plane C 



o r 3 



J(r) = f dzf(z) 

JCq 
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where 

z 3 

~ {z + z+){z + zJ){z - z+){z - z-){z + + —){z - - — ) ( ' A2 ' ) 



with 



1 



2 ± =( V ^+i v ^)-(l±VT^ 



r 

Here Co denotes the unit circle centered at the origin on C, \z\ = 1. There are four simple 
poles at z — Z-,~zZ, — z_ and — ~z~Z inside of the contour Co and the Cauchy residue theorem 



can be applied (see, for example, Chapter 4 in [37|] ) to obtain 

J(r) = 27ri{Res(/,z_) +Res(/,— ) + Res(/,-z_) + Res(/,-— )}, 
where we see 



Res(/,*_) = (z-z-)f{z) 

»-4 



2 7 ^pqVT^ (l-pr 2 )(l-qr 2 ) 



and Res(/, —-2-) = Res(/, Z-), Res(/, — z_) = Res(/, z_) = Res(/, We obtain 

t/ \ r 4 f 1 1 1 



24 yi^T 2 " [ 1 - pr 2 1 - gr^ 
The integral formula 



7 x arcsina . , . . 

ax- : = — , a < 1 (A3) 

o (l — a 2 x 2 )\/l — x 2 ay/1 — a 2 

is useful and we arrive at the result 

7T 2 

I = 7r(arcsin yfp + arcsin y/q) = —. 

It implies that fi p (v x ,v y ) given by (IT41) is well-normalized; f0^dv x J^dVy/ip^jVy) = I x 
2/tt 2 = 1. 

Similarly, we can also calculate the integrals 



/ / '' 2 

In. = dv T / an,,l 



OO /"00 II 1 

1* 



x = / GWx / ^^K/p+^Mi}/ , TTT7 TTT7 , 7T7 TT' 

2 

/oo /"oo -y^ 1 

^ y^ ^i MM/?<1} + + _ »^ + _ _ _ . 
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By the change of integral variables (1A1I) . we have 

I x = -2 2 ipJpq [ dr Jx ^\ L = 2 2 iqJpq [ dr 
Jo r jo 



with 



J x (r) = f dzf x (z), J y (r) = f dzf y (z), 



where f x (z) = (z + l/z) 2 f(z) and f y (z) = (z - l/z) 2 f(z) with llASI) . The Cauchy residue 
theorem gives 



4 ' 4 

ni r , . , m r 



2 2 (1 _ pr ) v ^372' 22 ( 1 _ gr ) v /i372- 

The integral formula (IA3[) and the fact Jq 1 c?rr/ y/1 — r 2 = 1 lead to the results 

I x = ^(arcsin^- y/pq), 
I y = 7r(arcsin y/q — y/pq). 

Since K x — I x x (2/7r 2 ) and K y = I y x (2/n 2 ), they give the expression ( |20l) . 

It is interesting to see that the above calculation of the integral I gives the following 
identity, 

I r 00 



dv x dv y /j, p (v x ,v y ) = drrfi(r;y/p) + drr/i(r; y/q), (A4) 

2 J — oo J—oo JO JO 



where /i(x;a) is the Konno density-function of one-dimensional quantum walk 12 . 
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VI - a 2 i 
fi{x; a) = — , =1 

7r(l — rjvfl 2 — x 2 



{\x\<\a\}- 
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